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ABSTRACT 


The work attempts to study and implement die sculptured 
surfaces on computer using various methods available in 
literature. The choice of parametric equation^ particularly 
B~spline system is justified for the surface modelling . 3D die 
surfaces have been generated by defining control points and by 
data fitting technique. Most of the required die sculptured 
surfaces are being generated using sweep technique of 2D curves. 
Open surfaces are modelled by sweeping 2D curves on two boundar 3 ^ 
curves to meet the requirement of forging or die cast surfaces. 
Extrusion die surfaces are obtained by cylindrical sweep 
technique considering the minimum requirement of the extrusion 
ratio to be constant. The joining of surface patches with and 
C continuity is implemented to improve the efficiency of surface 
desi gn. 

The system is based on the mathematical models 
suggested by various researchers engaged either in the field of 
surface modelling or in development of integrated CAD/CAM of 
dies. The software is implemented on HP“9000 work station using 
Star base graphics and C language. 
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CHAPTER 1 


INTRODUCTION 

In the 1970s^ the trend in advanced mass production 
technology was towards die based processing and moulding. 
Manuf acturing using dies and moulds offer numerous advantage over 
other manufacturing methods such as machining and joining 
processes. In medium to large production, moulding and forming 
methods become superior in terms of lower unit costs. In 
addition, material utilization is very high and rejection rates 
become very low. Dies are used in a variety of industries 
including pressing, forging, casting , rubber, glass, ceramics 
and plastics. Manufactured goods using cast parts range from 
airplanes to toasters and from automobiles to toys. 

Now a days, such industries as automobile, electrical 
goods and heavy machinery units use dies to manufacture a majority 
of their basic components. The highest level of technology is 
required for accurate and efficient die making, so that a wide 
variety of manufactured goods can be produced economically. 

The higher cost of overheads in addition to customer 
demand for variety of products at lower price and quick delivery 
make automation of die making industry an important priority. The 
spectacular growth of CAD/CAM - Computer Aided Design and Computer 
Aided Manufacture, which emerged within a last couple of decades 
as a result of revolutionary technology advancement of micro and 
digital computer, has helped to increase Engineering Productivity 



tremendously. 


There are inherent difficulties in designing extrusion 
and forging dies for making discrete parts. Traditionally, die 
design is considered more of an art than science. The analysis of 
irregular shaped forged dies or of streamlined extrusion dies for 
complex non-axis symmetric sections has been limited in the past 
due to complexity and high cost of die design and manufacture. 
However, the advances in the use of computers in Engineering 
Design and Manufacture have opened up new avenues in such 
application. 

In this context, we shall see later how the die surfaces 
can easily be represented on a computer screen before any analysis 
for modelling purpose which is done as the first step toward 
integrated CAD/CAM of dies. 



1.1 LITERATURE SURVEY 


Hajimu Kishi of Fujtitsu Fanuc LtdCll in 1970 defined 
and classified the sculptured surfaces as required for dies. He 
also suggested some technique to obtain die sculptured surfaces 
from sectional diagrams. Fanuc Ltd used this technique to develop 
a software FAPT Die-II for their CMC machines for design and 

manufacturing purposes . The software was designed for the desk 
top type CMC tape preparation system with graphic system 
considering the fact that expensive CAD/CAM systems are readily 
not adaptable to the economic need of an industry composed of 
small scale Job shops. 

Gunasekara and Hoshino C33 carried out many theoretical 
and experimental studies on extrusion dies having circular entry 
shape and regular polygonal exit shape and generated optimum die 
profiles using computers and used the data to produce engineering 
drawing and MC tapes for machining of EDM electrodes and 

subsequent electro discharge machining of dies. 

Dunkan [4] has reviewed the research in sheet metal 

forming using computer aided design for limit strain applicati ons. 
He has also described the application of this simple analysis to 
the problem of shape control in smoothly curved plane. 

Gunasekara et alC53 developed an efficient software for 
analysis of metal flow during forging of complex shapes using 

finite element method. They also developed a concept for design 



of streamlined dies of arbitrary product shape and manufactured 
dies for the extrusion of complex structures- Visualization and 
understanding of complex 3D geometry are enhanced by the use of 
shaded color graphics. 

GunasekaraCS] developed a software named STREAM which 
consists of three packages for modelling and design of shaped 
extrusion dies including re-entrant section and computer aided 
manufacture of the die surface. 

NagpalCS] used a quadratic polynomial function to design 
the streamlined extrusion die for complicated section. 

SabaroffC6] used computer technique for generating the 
spiral bevel tooth surface geometry by simulating gear cutting 
process and manufactured EDM electrodes for the forged dies. 

YangC7] proposes the conformal transformation for 
extrusion of helical shapes. However » the model offers limited 
extrusion shapes. 

Shev and LeeC7] developed an integrated CADXCAM/CAE 
method for optimum die surface design of 3D extrusion surfaces. 
Tension parameter » geometric coordinates and polynomial curves are 
used to define a bi parametric surface for the dies. The optimum 
die surface subsequently was manufactured by them using CAD/CAM 
software and NC machining. The proposed model claimed to obtain 
better extrusion performance. 



l.a SURFACE REPRESENTATION SCHEMES 


Though 'the representation and design of surfaces on 
computer is still a state of art subject ^ numerous ways are 
suggested for the correct representation. As such there is no 
universal form to define all real surfaces, hence the surface 
appropriate for a given problem depends on the application. Here 
are few examples of surface representation schemes carried on by 
different researchers . 

Polynomial cubic patches, Bezier patches, B-spline 
patches, Ferguson patches C83 are few surfaces patches capable of 
proving useful surface geometries. 

The Math Group at University of UtahElll devised a 
robust surface generation method applicable to arbitrarily located 
data using interpolation technique. 

Ver geest et al at Delft University of TechnologyC 123 
have suggested carrying out a software named SIPSURF for designing 
arbitrary surfaces and subsequent machining of these surfaces. 

Choi Cl 53 has described a method of constructing a 
smooth sculptured surface by using sweep technique. The software 
DUCT is partly based on this technique. 

Peigle Cl 63 in his survey ON NURBS has described how a 
rational B~spline can be used as a powerful tool for any type of 


surface generation. 



There are many sof*tware developed for surface modelling 
by various industries like Boeing, SDRC, Intergraph corporation 
etc These softwares cater to specific needs of the industry. The 
development of surface model is still under research and the best 
modeller coping with all the problems associated in a surface 
modeller is yet to come[153. 

1.3 OBJECTIVE OF THE PRESENT WORK 

The objective of the present work is limited to surface 
modelling for forge dies C open surface Dand extrusion diesC 
cylindrical surface D. The surface design on the computer screen 
shows the final shape of the product. After reviewing literature 
available in CAD/CAM of dies, it is observed that most of the work 
are meant for specific problem like extrusion dies or forging of 
helical gears etc with analysis of material flow and load required 
for the operation. 

In this work, an attempt has been made to model surfaces 
for dies as prescribed by KishiC23 which can be assumed to be the 
first step of CAD for dies. A complete CAD of die design would 
have been as shown in the flow chartCFigl . 12) . However, taking 
into account the volume of work involved in the analysis, it is 
emphasized in the modelling of various surfaces only and their 
presentation rather than detail analysis. Figl.3 briefly 

describes our work. 

We consider parametric curves and surfaces as two 
different tool to describe die sculptured surface. A basic 
surface can be formed either by supplying control points or by 
exact coordinates on the surface which could be fitted to a 



appr oxi mated surface. 


Sweep 'technique is used to describe most of the surfaces 

1 

as given by CS.3 . Similarly area mapping technique£73 is used 

along with sweep to define cylindrical surface for extrusion dies. 

An attempt is made to connect arbitrary shape surfaces 
with tangency or curvature continuity to facilitate the designer 
to describe blending etc. 

The software is implemented on the HP— 9000 work station 
using star base graphics and C language in unix environment. 




Fig 


1.1 FLOW CHART FOR DIE DESIGN AND MANUFACTURE 























CHAPTERS 


DIE SCULPTURED SURFACE AND MATHEMATICAL FORMS 

2.1 INTRODUCTION 

Die shaping can be divided into two general types: 2 
dimensional and 3 dimensional. These shapes can be characterized 
by planes, spheres, cones or cylinders etc so that NC milling 
machines or wire cut Electric Discharge machines can be used to 
machine these arithmetic shapes. 

Typically the die makers receive specification from the 
industry in the form of sectional diagrams for 3D sculptured 
surfaces. With these specifications the die makers are expected 
to form the intermediate sections and finished shape so that it is 
functional or pleasing to the eye. 


2.2 CLASSIFICATION AND DEFINITION OF DIE SCULPTURED SURFACE 

The die sculptured surfaces are defined by a combination 
of two kinds of section curves. One is called basic curve and the 
other is called drive curve. The basic curve is a shape element 
which governs the motion of drive curve. The drive curve is a 
shape element which moves with reference to the basic curve. Fig 
2.1 gives some example of die sculptured surfaces according to 
whether the condition in which the drive curve moves is parallel, 
radiate or normal[2]. Fig 2.2 shows sculptured surface for 2 
drive curves. Fig 2.3 shows some classification of surfaces 



depending on how the drive curve moves on the boundary curve. 


2.3 CYLINDRICAL SCULPTURED SURFACE 

Such shapes can be characterized by sectional curves on 
cylinders. The basic curve is Tixed and as drive curve moves 
along the basic curve the drive curve shape changes from 1 to 
2 Fig 2.4. The classification of cylindrical sculptured dies are 
shown in Fig 2.5 C23. The extrusion die surfaces fail into this 
category. 

2.4 MATHEMATICAL FORMS OF SCULPTURED SURFACES 

The geometry of a die may be approximated by a 
collection of plane faced polyhedra with straight edges. However » 
such a representation may involve hundred or thousands of faces or 
edges and would clearly be cumbersome to generate and modify. 
Hence it is necessary to represent the surface of a die as a 
collection of curved sculptured surface patches bounded by curved 
edges. However* irrespective of the method used to represent 
complex surfaces the geometry has to be modelled in mathematical 
-'form suitable for a computer. The following are some properties 
that mathematical model should haveClOl - 

1. Control points - may be on the curve or outside on a polygon. 

2. Multiple variable - a curve may have multiple values*!, e. for 
one particular coordinateCx* y or zD * there may be multiple value 
of other coordinates. 

3. Axis independence the shape should be independent of 

reference axis chosen. 



4r. Global or Local control — a Tine local adjustment to the surface 
should not affect the geometry as a whole. 

5. Versatility - mathematical models of curves should be versatile 

6. Orders of continuity - the designer may want to Join pieces of 
curves or surfaces maintaining the continuity of the geometry 
defined by the same control of points,!. e. the models should have 
well defined geometric properties. 

The above requirements have lead to the use of 
parametric or vector value methods for representation of curves 
and surfaces as given below. 

For curves 

PCu3 = C XCuD YCuD ZCuD ]. C2.1D 

For surfaces 

SCu,vD = [ XCu,vD YCu,vD ZCu,vD3. C2.2D 

2.4 SURFACE PATCHES 

Computer aided design objects are not solely concerned 
with pictorial properties of shape but with providing a geometric 
model or geometric data base which can then be used for 
design, analysis and manufacture. Some well-known surface patches 
used in various industries are bilinear surface, Coon's surface. 
Lofted surface, Bezier surface, B-spline surface and rational 
B-spline surface patchesC93. 

However, it is not easy for a designer to to use Coon's 
surface patches or Lofted surface patches to define any arbitrary 
shape. Bezier or B-spline systems are widely used for describing 
surfaces of automobile, ship hull, aeroplane body etc. In this 



work w© have used B sli.ne or Bezier curves for defining surf ace 
patches. 

The parametric equation of bezier surface is - 
n m 

PCu.vD = £ E P. .B. CuD B. CvD CB.SD 

i=o j=o 

where the Bernstein Polynomials B's are given by 


B CuD = '^C, u'- Cl-uD 

i.,r. 

B, CvD = Cl-vD 

j 


Tr»-j 


ca. 4D 

ca, 53 


wi th 


t ! < ri"“ t > ! 


and 


ca. 63 


I j ! < m- j >! 


^ are the vertices of the defining polygon net. n and m 
are one less than the total number of polygon vertices in the u 
and V direction respectively. u and v are two parameters. 


The parametric equation for B~spline is ~ 

n m 


SCu.v3 = Y. 


Y P. . N. Cu3 M Cv3 

. ^ ^ Wk I 

t=o j=o 


ca. 73 


where N Cu3 Cv3 are B-spline basis functions in 

the bi parametric u and v directions respectively. 


N. = I 


if X < u < X 

t t. + 1 

else otherwise. 


C£. 8D 


Cu- X.D N. . CuD 
Ni.k = ; — — + — iii 


CX - N, , Cu:> 

\.+±,k 


ca. 9D 


X. , - X, 

v+k-l V 


X. , - X, 

1, + k t+l 


and 


M 


= r 

I O 


if y.< V < y. 

1. t +1 


else otherwise 


Ca.l03 




where P, . are the vertices of defining polygon patches. For 
quadrilateral surface patches the defining polygon net must be 
topologically rectangular. The Indices n and m are one less than 
total number of polygon points in u and v direction respectively, 
k and 1 are order of the curve in u and v direction. In the above 
equation 0/0 is taken as O. 

X, and Y. are called knot vectors which are important in 

X. h 

deciding curve and surface shape. There are 3 types of knot 
vector namely - Open uniform^ Close uniform and Mon~uniform[83 . 

When open surface is required open uniform knot vector 
is used and given by ~ 

X.= O i < k . 

t 

X.= i-k+1 k< i < n , 

i 

X.= n~k-^2 i > n . 

i 

The parameter u lies in the range of [0>n~k+23. 

When a closed surface is required the close knot vector 
is used and given by - 

X. = i . C2.153 

I 

and Ck-2D polygon vertices have to be repeated at end/or beginning 
of each curve. The parameter lies in the range of Ck~iD to n~k+l. 

Nonuniform knot vector is useful in connecting separate 




segments of a curve and defining quadratic curves like circle, 
ellipse etc and are decided by various mathematical forms[8]. 

Fig 2.6 and Fig 2.7 are two examples of open and close 
B-spline surfaces respectively. 


NONUNIFORM RATIONAL B-SPLINE: 

The rational form of B— spline offers a common 
mathematical equation for representing and designing both standard 
analytic shapesC conics, quadratics, surf ace of revolution etcD and 
free form surfaces precisely. The parametric model is given by- 


n rri 

E E V. . P, . N. Cu3 M Cv:) 
. " '-'J '"J j'l- 

v = 0j=0 

SCu,vD = 

n m 

EE W . N. CuD M. CvD 

, ^ Wj t,k 


C£. 16:) 


All notations are 
W is a weight function 
vertices. When W is 1 for 
ordinary B~splineC83. 


the same as given for a B-spline surface, 
associated with each polygon control 
all i,J rational B-spline reduces to 


Z.S WHY B-SPLINE IS BEST FOR DIE SURFACE ? 

In the present work, B-spline is used as design curve or 
surface because it shows all necessary properties mentioned 
before. It has advantage over Bezier system as it allows local 
control which means the influence of any single vertex of polygon 
net is limited to only to -k/2 span on each curve. Again 

when the number of defining polygon net vertices is equal to the 



order of each curve B-spline curve reduces to Bezier curve. 

In addition^ it has some more advantages given belowClS]- 

1. It can be manipulated in number a of ways such as by tweaking 
polygon net vertices, order of the curve, using multiple polygon 
vertices , type of knots , multiple knots etc . 

2. Evaluation is reasonably fast and computational ly stable. 

3. It has; a powerful geometric tool kit;i.e. scaling, rotation, 
translation, parallel or perspective transformation can be carried 
out on the geometric database. 

4. B-splines have clear geometric interpretations, making them 
particularly useful for designers who have a good knowledge of 
geometry especially descriptive geometry. 

5. The B-spline surface algorithms is mostly supported by work 
stations. 

6. The non-uniform rational form of B-spli neCKURBD are genuine 
generalizations of nonrational and rational B-spline system as 
well as Bezier curves. But it needs more storage and demands 
careful choice of weights associated with each polygon vertices. 
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TyfM 0f icwlptund iurfm 


1st digit m No, of be, 2nd digit » No, of drive cur 


condition in which dc moves on be 
Typical Shapa of Sculptyrad Surface 


Sculptured Surface 205 
(connection) 

Sculptured Surface 1 1 1 
(panllel) 

Sculptured Surface 112 
(radiate) 

Sculptured Surface 113 
(normal) 

Sculptured Surface 121 
(panllel) 

Sculptured Surface 122 
(radiate) 

Sculptured Surface 123 
(normal) 

Sculptured Surface 21 1 
(parallel) 

Sculptured Surface 212 
(radiate) 

Sculptured Surface 213 
(normal) 

Sculptured Surface 214 
(perpendicular) 

Sculptured Surface 221 
(parallel) 

Sculptured Surface 222 
(radiate) 

SculptHied Surface 223 
(normal) 

Sculptured Surface 224 
(perpendicular) 



Fig, 2,3 Classification of Sculptured Surface'^ 
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Fig. 


E. 4 Cylindrical Die Sculptured Surface 
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Fig. 2.5 Classification of Cylindrical Sculptured Surface 
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Fig. 2*8 Exan^le of Close B— spline Surface 


Fig. E. 7 Example of Open B-spline Surface 
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CHAPTER 3 


SURFACE FITTING FOR SHAPE DESIGN 


3.1 GENERAL INTRODUCTION 

Sijrjra.ct& fi. 'L'Li. rxQ t/©chnic}Ui^s Irsive bosn d^v^iopoci in 3 . wid© 
variety and for a great number of different applications through 
interpolation or approximation scheme. The fitting methods 
usually vary with respect to the constraints which may be imposed 
such as boundary conditions^ smoothness criteria or local 
deviation tol erance. 

The interactive surface design can be carried out on 
Bezier /B*-spline system doing operations like control points 
tweaking etc as described in Chapter 2. However, the designer has 
limited control over the actual surface by changing control 
polygon vertices, hence editing is not suitable for faring curve 
or surfaces. Again for B-spline surfaces, the tensor product 
concept hampers control point addition and deletion because a knot 
is associated with a row of de Boor control pol ygonC Basi s 
function!) [16], 


Here a fitting of data point scheme is described which 
could be very useful for open or close die surface development. 
But the scheme is based on global method;!, e. change of any data 
point affect the surface as a whole however small may be. As a 
result the modelling is very fast for interactive process as long 
as the number of data is not excessively large. 

The fitting method is implemented using B-spline systeffii 
and the surface is optimized with respect to deviation from data 

23 


points in a least square sense which means the surface will 

deviate from any data point by less than some user specified 
distance. 

3.2 B-SPLINE SYSTEM 

The parametric form of B-spline surface is given in 
Chapter 2 EqC2. 8D to EqC2. 123. The equation can be written by 3 
coordinate functions in two common real variablesC parameters 3 u 

and V as 


nu nv 


XCu3 

= E E 

t =Aj=l 

Px 

N. , Cu3 

N C vD 

J^kv 

C3. 13 


nu nv 





YCu3 

= E E 

t =ij=i 


H, C uD 

i.,ku 

N C v3 

C3. 23 


nu nv 





ZCu3 

= E E 

i = ij= 1 

Px. . 

N, , Cu3 

i.,ku 

N., Cv3 

j^kv 

C3. 33 

and P, . = 

L. i 

C Px, 

Py. . 

Pz, . D 


C3. 43 


<■.0 WJ L,J V,J 


are three dimensional control points in Cnuxnv> size matrix, nu 

and nv are total number of control points. N and N are 

i/ku j,lcv 

B-spline basis function given by EqC2.93 and EqC2. 123. ku and kv 
are orders of the curve in u and v parameter direction. 

The possible surface topology are : 

1. Surface open in u and v direction C plane topology3. 

2. Surface open in u direction and close in v direction or vice 
versaC cylindrical topology3. 

3. Surface close in both u and v direction C toroidal topology3. 

The knot vector is used depending on whether surface is 
closed or open in that direction as given in EqC2. 143 or EqC2. 153. 
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SURFACE FITTING TECHNIQUE 
3-3 MINIMUM REQUIREMENT 

1 . The data to be fitted must be ordered in some way, such 
that contours or other characters tic curves follow on or near the 
intended shape. 

c-.If data originates from some measuring system then 
samples are taken in series or some gridding pattern. Normally any 
industrial design item does not have random coordinates. 

3.4 PRINCIPLE 

The surface fitting technique is done by repeated curve 
fitting. Suppose the data points are arranged as M series of N 
rit-s , then each series can be approximated by a curve 

represented by K coefficients. Now there are K series of M points 
available, each of which can be approximated by a curve 

represented by L coefficients using exactly the same fitting 
procedure. The resulting CKXLD coefficients are exactly the 
control points of the intended approximate surface. 

3.5 DESCRIPTION 

The following assumptions are made for the input- 
1 . The data consists of 3D cartesian pointsC®g x y z coordinates!) . 

S. The data points must be grouped in M series. 

3. Each series may consist of different number of data pointsCNl ; 
i = 1 , . . , Mi) . 

4. The order in which data points occur in a series are 
relevant, as is the ordering of series themselves. 

On the basis of these data a B-spline surface of order 
ku,kv will be found that approximates each point IX.j within 
tolerance T= CTx Ty TaD in a sense that for each D. ^ ^ there exits a 
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parametric value pairCu' v’D such that 

ISCu’.v'!) - Dx. I < Tx 
»• . J 

|SCU' .v*:) - Dy . I < Ty 

^ ^ J 

|SCu’ , v'D - Dz, J < Tz 

1 r J ' 

Tx Ty and Tz may be assigned some positive real value. T=CO O OD 
would result in interpolating B-spline surface that passes through 
all data points. 

3.6 CURVE FIT 

Initially the x-, y- and z- components of D are 

v--j 

fitted independently by single valued B-spline function XCu3, YCuD 
and ZCuD of orders ku defined on a knot vector as given in 
EqCE.143 or EqCB. IS!). CNote: N> ku or 1< N < ku, where N is total 
number of data is not pre requisite for the fit. These 

under deter mined cases will be solved appropriately as described 
later in this chapter^. 

The B-spline function XCuD with the lowest possible 

number of control points is fitted to the data. Ini ti ally number of 

control points nu is assumed to be ku.If XCuD fails to approximate 

the value Dx within tolerance Of 1/2 Tx then nu is increased by 1 
i 

and fit is retried. So when nu = N, exact data interpolation is 
obtained with all data points are on the curve. Thus a 

satisfactory fit is obtained with smallest possible nu existing in 
the range ku S nu < N. 

Fitting of 

nu 

XCu3 = r Px. N. , CuD 03.5!) 

i s=l 

to the values Dx. , . . . ,Dx proceeds as follows - 

i. n 
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The knot vector for u parameter Is setup for B-splipe curve 
of order ku and of control points nu. To each data point .a 

parameter value ti Is assigned in the range £ umir, Umax! , -Ih. ^ fit 

function is formulated of a function XCul to data points Dx such 

j 

that 

nu 2 

E W.CXCt.3 - DxD 


C3. 6D 


is mi ni mal . 

Wj is a weight factor associated with each data point Dj. 

But, the parameter value tj has to be found out in some 
way such that the parameter is distributed in some average pattern 

within the range [umin.umox]. 


3*7 CHOICE OF t 

j 

A minimal requirement for tj is that tj+i > tj for 1< j 
V N and there is at least one t in each nonzero interval 
span C ^ 3 for i= ku. . »nu. This still leaves many different 
possible choice for tj. 

Here, we consider distribution pattern suggested by 
Ver geest C 123 which gives regular spacing of t values within the 
parameter range. Suppose that a vector for u parameter has been 
set up for order ku and nu number of control points then there are 
nu-ku+1 knot intervals. 

First, nu locations s.are determined from- 


U, + u u 

t-t-2 f+*k-l 

^5 , I, „ 'i i.T ' jji ,1 1 1 , ' r ir n ^ JL {ji , ■ , ITlU 

1 k u ^1 


The s are scaled and shifted to produce 

i 

will be at extremes of u. and u 

ku nu-t-i 

s " = s /r - d , 1= 1 , . . nu . 

i t 


C3.73 ” 

s* and s^ which 
1 nu 


C3, 8D 
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where 


, _ __ C3. 9D 

nu+i ^ku 

s 

d = ~ C3.10D 

To'ta.l number of daba. poin“ts N is larger “than or egual *to 

nu ; therefore if 
Cn ~ 13 

A = lb C3.11D 

then, 0 ^ A < 1 . 

Finally, tj is taken as follows: 

tj = s’ + C AC J -1 3 +1 -nDCs ’ -s’D C3. 13D 

^ m iD+i m 

where 

m = INTCACJ-1D> + 1 C3. 14D 

and INTCaD is the largest integer less than or equal to 
a C floor function of C -language^ . 

The above distribution guarantees that ti and tn are 

located at the interval bounds and there at least lies a nonzero 
value within each interval .which is supposedly an appropriate 
condition . 

If the weights EqC3.6D are taken identical and 

equal to 1, then with the knot vector and tj parameter described 
above it can be verified that there exists an unique set of 

Px,..Px that solves EqC3. 6D. The vector form of EqC3. 6D can be 

nu 

written as - 


N , CtiD 

1 , ku 

,kuCt23 
N . CtND 

t , ku 


N ,CtiD 

nu » K u 

N , CtzD 
nu , ku 


N , CtN!) 
nu , ku 



E)x 

1 

Dx 

2 

Dx 

N 


C3. 14D 
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sum of squares of 


The expression 3.14 has minimal lengthCi.e. 
components has the minimum valueD. The above matrix is of NXl size 
and has many zeroes if ku << nu << N due to local support of 
b-spline function. 

^ C3. 15D 

The above matrix can be solved as given below. 

Differentiating EgC3. 6D for least square error and 
equating to zero gives the solution as 

[b’’]CB][P3 = Cb'^]CD] C3.16D 

where 


CB3 is the basis function matrix given in EqC3.14D. 

EqC 3 . 1 6D becomes 

CA3CP3 = Ce] C3.17D 

How, CA3 is a positive, semi -def ini te, symmetric matrix. 
So this makes the system suited for solution by Cholesky 

factorisation methodC183. The solution gives the control points 

Pxi, 

The same procedure is followed to obtain yCu!> and zCuD 

to approximate Dy and Dz within tolerance i/zTy and i/zTz 

respectively. But if one or more components deviate too much from 

the data ,nu is incremented for all 3 fits with XCu3 , yCuD and 

zCuD. At this point we have obtained a set of control points PCiD 

to draw a B-spline curve of order ku. 

nu 

PCuD ^ T Pi N, , CuD C3. 183 

t , k u 

i i 
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3.8 RE-REPRESENTATION OF B-SFLINE CURVE 

For each of the M rows of data points a b_spline curve 

r ,C uD can be fitted in the way described above. If different nu*s 
^ j 

are obtained for different rows of data, a choice has to be made 
for one single value to keep knot vectors identical for all rows. 
The following choice is made here - 

numax = Maximum of nui,nu 2 ,. . . . ,nun^. C3. IID 

3.9 UNDER-DETERMINED SYSTEM 

1. If numax > Nj the system gets under -deter mi nedC 133 . The 

following method to represent the curve removes the above problem. 

numax locations e. . on curve r C uD are evaluated at 

j 

parameter value ti, i= i. .numax and ti are newly determined for 
numax control points as per EqC3.73 to EqC3.13D. Then locations 
e . are regarded as numax data points to B-spline curve of order 

ij 

ku to be fitted or interpolated. 

The above method keeps the original fitted curve 
principally unchanged except for some extreme cases with large 
di sconti nui ti es . 

2. If N. < ku the system gets under deter mined - 

To handle the above problem , the order is increased to Hand 

nu is initialized as N,. Then the whole procedure is repeated as 

J 

described above to find numax control points of the curve. 

3.10 SURFACE FITTING FROM FITTED CURVES 

We have now obtained M B-spline curves r.CuD of order ku 
with numax control points E_ Ci=i. .. numax, j= i. ..mD. How we 
regard these control points as numax rows each consisting of M 
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points. For each row i we determine B-spline curve r CvD of order 

i 

kv that approaches M points E,, within tolerance 1/2 T and with 
minimum number of points. Then the same procedure as for u 
parameter is followed to find nvmax control points for B-spline 

surface b.<-U,vD of order ku »kv. As the curve r CuD deviates i/2T 

j 

from original data points D while curves r CvD are closer than 

i 

1 / 2 T t-o cont/Tol poi nt^s of curve r^CvZ)> "the surface approxima.'te t.o 
original data by tolerance of T C 13 ] except for perhaps some 
extreme cases. In the current i mpl emententati on there is no 
explicit distance check between D and S 

3*11 SURFACE FITTING FOR DIES 

In the past, the die surface model was based on hit and 
trial method. Despite the computer provides excellent facilities 
for designing surfaces, in most cases computer graphics proves 
uneconomical for dies because of the effort involved in developing 
programs in describing the geometry of the shape. Although the 
parametric surface can easily be made by means of CAD technique, 
it is obvious that mathematical expression can not exactly match 
the real sculptured surface. In order to get more adequate 
parametric surface, more detail data points should lie within each 
patch and the surface should be fitted to those points. 

The surface fitting technique described above may be 
useful for following two purposes. 

1 . RKFUOfttlCTION OF WORN OUT ME SURFACE 

When A die surface is worn out. the surface may be required 
to be reproduced with some alteration. The old surface can be 
measured by Coordinate Measuring Machine by moving its probe on 


31 



the surface; for example 100 data are measured in 10 rows. Taking 
probe error into account, the original surface can be fitted by an 

approximated surface for further manipulation and subsequent 
manufacture. 

2. PRODUCTION OF A NEW DIE SURFACE 

In spite of development of CAD/CAM system the designer 
still needs models made of wood, clay, plastic or leaflets because 
their design idea of product shapes is difficult to define in 
geometric expression and visualize on a 2D picture plane. 

After making those models by artists, the surface is to 
be generated for NC milling or other computer controlled 
production processes. The surface fitting technique may help for 
that in a big way. 

The fitting technique is particularly useful for defining 
dies for forging or pressing car body, complicated mould 
surfaces, cores in die casting etc. 

3.11 I MPLEMENT ATI ON 

The fitting of surface has been implemented in HP work 
station. The data can be simple fed manually key board or through 
a text file. To design a die sculptured surface, the data can be 
obtained by a CMM -coordinate measuring machine from a clay model 
or by 2-D or 3-D digitizer and arranged in a grid pattern for the 
input . 

3.12 CONCLUSION 

Because of the B-spline basis function support, the 
above least square fitting technique is much faster than least 
square solver E04FDF of Nag library which uses general least 
square method and therefore our method is more suitable for 
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interactive surface development. However , the control points found 
t’y are not easily interpretable generally. Smoothness and 

closeness to the data is controlled by simple tolerance values 
presently and interactive manipulation of data point is not 
implemented in our work. Nevertheless it is expected that above 
insthod of surface fitting to be of immense help for die surface 
model ling. 
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CHAPTER 4 


SWEEP SURFACE MODELLING FOR DIES 

4.1 GENERAL INTRODUCTION 

A smooth sculptured surface that is best described as 
trajectory of cross-section curves swept along profile curves is 
called a sweep surface. It is easy to describe aesthetically 
appealling sweep surfaces by specifying only 2-D cross 
sections. The surface SCu.v:) can be represented as a series of 
simple transformation and blending of cross section curves and 
profile curves. In the proposed model , a sweep surface is 
regarded as a tuple consisting of cross section curves, profile 
curves and a sweeping rule. The sweeping rule considered are 
parallel sweeping .rotational sweeping, spiral sweeping and 
synchronized sweeping. 

Using the above rules, most of the sculptured surfaces 
as prescribed by KishiCS] can be described. 

4.E DESCRIPTION OF SWEEP SURFACE 

In engineering drawings, a sweep surface is 
usually specified by its profile curves and a number of cross 
section views as shown in Fig. 4.1 . Presented in Fig. 4.2 is 

the 3D view of the sweep surface. The base plane corresponds to 
Fig 4.1Cb:) and the first left cross section plane represents A-A” 
section in Fig. 4. ICdil.The section plane can be assumed a solid 
plate on which section curve is drawn and the plate is moved 

in space so that b^CvD passes through p^ while the plate is 
vertical . 
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n a unique 


A sweeping rule has to be specified to obtai 
motion of section plane. If the plate moves in such a way that 
the flow rates on the two boundary curves are the same, e. g. 
section plane intersect at the same percentage of the length of 
the boundary curves, then the surface is called synchronized 
sweep. 

The above sweeping rule guarantees that first boundary 
curve b^Cv3 always passes through p^ of the moving plate, but the 
second boundary curve b^CvD does not in general intersect the 
plate at p^as depicted in Fig 4.3. Thus it is necessary to 
correct the original section so that the initial p is brought to 
the current intersection point. The corrected trajectory of the 
first section curve defines a sweep surface. By repeatinq the 
same steps for second cross section curve another sweep surface 
is obtained. Finally the two sweep surfaces are blended 

according to a suitable blending rule to obtain the final sweep 
surface. 

4.3 MATHEMATICAL REPRESENTATION 

A mathematical model of surface is a mapping from a 3D 
domainCof parameter u,v3 to an space. Our goal is to obtain a 
mapping function rCu,vj from the description of two boundary 
curves and two section curves . The following steps give the 
mathematical replica of the descriptive procedure of the previous 
section. 
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4.4 COORDINATE FRAMES AND CURVE EQUATION 


The two bounded curves are represented as parametric 
curve equations with respect to a "profile parameter" v: 

bCvD=CxCv3)yCvD2CvDD o<v<-i 

b^CvD = C x^CvD yCv3 zCvD , o < v < i C4.2D 

The requirement here is that mapping from v e[o. i] to bCvZ) 
is regular and non-singular, that is one to one mapping is 

requi red. 


Let. I VC vD in Fiq. 4.3 denoLe Lhe vector from b CvD Lo 

1 

X— y plane. The angle ©CvD between x axis and 2D vector 
IVCvD becomes function of v and is given by, 

eCvD = a tan2C yCvD ~ yCvD , xCvD ~ xCvD D C4.33 

o 1 o 1 

In order to have a well defined intersection angle eCvD 
the projected images of two boundary curves are not allowed to 
intersect with each other. 

Each section curve is represented as parametric curve 
sCuD with o < u < ± Csection curve parameter u3 . Also specified 
in fig 4.5 are coordinates of guide point G through which 
passes. The intersection vector IVCvD appears on the section 
plane and the angle (j) between x~axis and the intersection vector. 
The guide point is as follows- 

G = Cq g 03 = sCo3 C4. 43 

But it is always possible to define the section 
coordinate frame such that G = CO O 03 and cj> = ISO'^which makes 
the transformation matrix much simpler. 


4-5 SWEEP TRANSFORMATION 

The following sequence of coordinate transformations 
bring the section plane of x-y frame to a 3 dimensional location 
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of a section curve[14]. 


1. Translate so that guide point moves to origin. 

E. Rotate - through an angle <|) around z axis. 

3. Rotate - through an angle -90 degree around x-axis. 

4r. Rot-a.!.© fcjy oCvD around z axis. 

5. Translate so that G moves back to b CvD 

o 

The 3 rotations from the core of the sweep 
transformation can be conveniently expressed in matrix form. 
SweepC4>.eCvDD = RotCz.-cjO RotCx, -90D RotCz. eCvD!) 
which can be expressed in mathematical form [8] 


cosc|> -sincfc O 
sincj) COSO O 
O O 1 



coseCvD sineCvD 
-si neC vI) coseC vD 
_ 0 O 



C4. 5D 


cosc|> coseCvD 
= sincf) coseCvD 
-si n©C vD 

Where c|> and eCvD are intersection angles on the 
section plane and base plane, respectively. 

If <|> is chosen as 180*^, the sweep transformation 
matrix simplifies to- 


C0S(J> sineCvD sin<j) 
sin(|> sineCvD — cos(j) 


coseC vD 


C4. 6:) 


SweepC eC v3 3 


-coseCv3 -sineCv3 
O O 

jsin©Cv3 cos©Cv3 



C4. 73 


In other words the sweep surface PCu,v3 is obtained by 
sweeping a section curve sCu3 along a boundary curve b^Cv3 and is 
expressed as the following sweep transformation. 

PCu,v3 = CsCu3- G3 SweepC 4>,eCv3 3 +b^Cv3; O < u,v < 1 C4.83 

If guide point G becomes the origin and intersection 
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vector lYCvD is in minus x-axis direction Cj;e.G=CO O OD and c[) = 
100 I) » then sweep surtac^ ecjuation becomes 

PCu,vD = sCuDSweepC©Cv3D + b^CvD ; 0 < u.v < 1. C4.9:> 

4.6 CORRECTION FOR SURFACE BOUNDARIES 

As described earlier , the end of section curve in 

general is away from boundary curve b^CvD, the curve has to be 

brought to b^CvD by some way for the completion of the surface. 

The following method is suggested by Choi and Lee[14]. 

First, the section curve is scaled so that distance 

between the 'start point' pCo.vD and 'end point’ pCi.vD becomes 

egual to the distance bistween two boundary curves. Second, a 

section curve is obtained so that pCo,vD coincides with b CvD. 

o 

Third, another section' curve is obtained by moving PCi.vD to 

b^C vZ) ; and finally blending two two resulting section curves by 

using a blending function. The scaling factor fCvD is defined as 
I b C vD - b C vD I 

fCvD = i C4.10:! 

I sCoD - sCiD I 

The corrected intermediate section curve qCu,vD is then 
expressed as a blend of two scaled and translated intermediate 
section curves ; 

qC u , vD = ocC uD< f C vD [ pC u , vD - pC O , vD 3 + b^C vD > 

+ BCuiXfCvDCpCu, vD - pCl , vD 3 +b^C vD > . C4.11D 

where oc and B are two-bldnding functions. 

4.7 BLENDING 

Similarly, two' section curves can be blended by a 
suitable blending function. The equation of the blended sweep 
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surface is expressed as : 


rCu.v:) = ocCvDq^Cu,vD + fiCvD q^Cu,vD 

= MCvDf^CvDCs^Cu:) - otCu Ds^CoD - BCuDs^C iD ] SweepC^)^ , eC vD D 
+ eCvDf^CvDCs^CuD - ocC uDs^CoD - BCuDs^CO] SweepC 4>^ , eC vD D 
vD BCuDb^Cv])!! ^ U, V < 1 . C4. 120 

where s^CuD are equations of section curves i = o,i ; fCvD is 
the correction factor as given in eq. 4.10 . c|> ,SweepC<j), eCvOO 

are as defined earlier. aC*0 and BC*0 are two blending function 
chosen such that ocCo) = BCIO = 1, odClO = BCOO = O, and <xCuO+BCuO 
= txC vi) +BC uD = 1 . 


The choice of the blending function otC*0 and BC*0 in 
eqC4.1SD affects the quality of the surface. The simplest choice 
is linear blending function: 

ocCuD = 1- u ; BCuD = u; C4. 13:> 

Another blending function can be chosen called Her mite 
blending functionClbl ; given as ; 

ocCuD =1 - 3u^ + 2u^ ; BC uD = 1 - ocC u^) ; C4. 14Z> 

The above blending function is suitable for providing 
tangent continuity and therefore can be used when a sweep surface 
is constructed from network of boundary and section curves. 

The above model describes a synchronized sweeping 

surf ace. 


4.8 GENERAL SWEEP SURFACE 

For the general sweep surfaces, the concept of guide 
curve is introduced. A guide curve could either be spine curve 
or boundary curve. The boundary curve should not necessarily lie 
on the surface. For example, if the section curve is a closed 
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curv* .then the centre or cross section of the closed curve cen 
be teken as guide curve. A parametric form of guide curve is: 

gCvD = CXgCvD YgCvD ZgCvD3 
Derivative of the guide curve is ; 

dgCvD/dv = CXgCVD YgCVD 2gCV:)3 ^-4 

The intersection angle eCvD under different sweeping rules are 
gi ven as foil ows ; 

1 . ©CvD - eo C Par all el sweepD C 4 . 17D 

c. . OC vl) — atandCYgCvD ■" ay » XgCvD “ ax 3 

C Rotational sweepD 

a = tax ay oD ; coordinate of the axis of roatation. 


3. ©tvD = a tanSCXgCvD .-YgCvDD 


C Spi ned sweepD 


C4. 19D 


4. eCvD = a tanaCYoCvD - YiCv3, XoCv]) - XiCv53 C4. 20:3 

C synchronized sweepD 

For a given value of v the profile parameter and guide curve 
gCv3 the sweep transformation of Eq.C4.9D can be modified as: 

PCu,v3 = CsCuD - GDSweepC(|> .eCvDD + gCvD C4.21D 


4.9 GEMERAL CORRECTION TRANSFORMATION 

The points on the intermediate section plane at which 

the two boundary curves intersect with it are determined by 

numerically solving the following equationsC 18] . 

b C s]> . nC vD = qCvD.nCvD C4.22D 

o ^ 

"b” C tD . nt v3 ="^vD.'nCvD C4. 233 

o 

where nCv3 is the normal vector of section plane. 

EXAMPL.E OF SOLUTION: 

The parametric form of B-spline or Bezier curve are given by 
Cohen and Riesenf eldt 163 . For example Bezier curve can be 
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deduced in the form - 

pCuD = [U][N][G] 

wher e , 


C4. 24rD 


r I n r 

t U3 = t U U 


Nl+i ► j + 1 


u iJ 


C4. 25D 


J 


, IC'CjD C^'^Cn-i-j D- 

I O < i+j < n 

o I O otherwise. 


[G] = C p p 
o 1 


P ] 

n 


C4. 26D 

C4. 27D 

Usi nci ©‘qua. ‘Lions EqC4.21D "to EqC4.27D^ a polynomial function 
can be formed whose roots may be real and distinct, real or equal 
or complex root£>, S^o, the solution can be found out by Oraff *s 
root squaring method[18]. 

Let s and t be the solutions of the equations. Then 
scaling factor fCvl) for general sweeping rule is given by 

fCvD = |b Cs*:) ~ bCt*D|/| sCo2) - sCO | C4.£S:> 

The final generalized sweep surface can be 
mathematically expressed in the form: 

r C u , vD cx< vD f C vD [ s C uD cxC uD s C u3 ~ BC uD s Cl D 3 SweepC d) , OC vD D 

o o o o 

■+-BC vD f vD [ s^C u3 - cxCuDs^CoD - BCuDs^CiD 3 SweepC <|>^ .eCvDZ) 

+ [cxCuDb Cs'^D + BCuD b Ct'^D3 ; 0 < u,v < 1. C 4.393 

o 1 

In summary, a sweep surface has one guide curve, up two 
boundary curves and one or two sweep curves. 

For completeness, it may be noted that the die sculptured 
surface of fig 2,. 2 can be generated using sweep technique and 
different combination of boundary and section curve or curves. 

The above surface model technique could be used for 
forging, moulding and casting dies. 
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without twisting or bending. 

The ba.sic concept of extrusion die geometry is shown in 
figC4. D. The sector of the billet has been mapped onto 
corresponding sector on each side with same extrusion ratio being 
maintained. Hence, 

Area OPQ / Area OP'Q’ = R = Area OPR / Area OP’R' 

where R is the extrusion r ati o. GunasekaraC 33 suggested that spline of 
any order can be fitted to the entry and exit sections to obtain a 
streamline dies. 

But, the area mapping technique fails for a re-entrant 
type die geometry as shown in fig.C4.8 D. Area mapping technique 
is transformed to perimeter technique by use of Stokes 
theoremClO] which means the perimeter of the billet has to be 
mapped to the perimeter of the extruded product. This solution 
is more generalized and applicable even to re-entrant section. 

In our work the following method of cylindrical sweep 
surface has been used to describe extrusion die geometries: 

A. SYMMETRIC CASE: 


4.112 Product Section : 

The coordinates of extruded product is decided 
first in two dimensional Cx-yD plane. If the product section has 
all straight line edges then a B-spline of order 2 is fitted to 
the control points to form a close curve. If the section has 
curves a higher C3rd or 4th3) order spline is fitted. 

4.113 Symmetric Axis: 

On the product boundary axis of symmetric has to be 
decided. The axis of symmetry is taken to be the x and y axes 
through the centroid of the product. Centroid is decided as. 
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Xc 


1 iXl + I 2 Xz + ... +lnXri 


C4. 30D 


1 1 + 1 2 + . . . . +1 n 

1 lYi + I 2 Y2+. . ■ +lnXn 
1 1 +1 2 + . . . + 1 n 


C4. 31D 


where li,l2 . . etc are the segment length on the product section, 
and Xi.Yi . .etc are the coordinates of the centroid of each 
segment measured from the centre. 

4.114 Die Parameter 

The die parameter such as billet diameter and die length are 
decided from analysis for optimum design and can be fed as input. 
The billet diameter is drawn taking centroid of the product as 
the centre. 

4.115 Surface Construction 


1. An adequate number of char actersti c points on the exit curve 


are selected. 

S. The mapping points on the billet are determined by connecting 
a straight line between the centre point CXc^YcD and the 
character isti c points on the exit boundary. 

3. V-curve construction C Axial di recti onD 

It is assumed that the axial profile of the die are 
streamlines which can be designed with flexibility. To draw a 
streamline, two more control points are decided between product 
section's characterstic points 'Pi' and corresponding mapped 
point 'Qi.* on 'billet di ameter ' . the derivative of the B-spline at 

Pi and Qi are 


p[ = 3CPi - EiD . Qi’= 3CBv - Oil 

where Bi and Ei are the Z coordinates of 

corresponding to Qi and Pi CX.YI) coordinate 


C4. 321) 

the control points 
figC4.9D. Assuming 
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axial entry and exit .Ctangentis zeroD z coordinates are taken 
as 1 /AC die length!) from Pi and Qi. The streamline constructed by 
Qi,Bv,Ei,Pi. are shown in figC4.9 

Now, four control points are obtained for each stream line. 
The stream line can be swept on two cylindridal curves for 360° 
to produce a extrusion die geometry. This can be done in simpler 
manner. The four points are decided through all characterstic 
points and a close B-spline surface can be fitted through these 
points. 

4.116 Unsymmetrical Case: 

For unsymmetric section with re-entrant shape the perimeter 
may not map exactly as described before. We suggest for those 
cases that node 1 of the product cross section should start from 
a point on the billet which intersects the positive y-axis and 
the rest of the procedure follows as above. 

However, a detailed analysis has to be done on how the 
material flow occurs in these types of surfaces to decide the 
smooth stream lines. 

Some of the examples are given at the end for illustration. 
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c c 

Fig.A.l Description of sweep surface‘s 



Fig. 4. 2 3D view of sweep surface'"' 
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Fig. 4. 5 Section coordinate system for 
synchronised sweeping. 
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)ng.4.6 Swacpirig ruUs; (o) porallcl Fig. 4. 7 Section coordinate 

(b) rotational (c) spined (d) synchronized.'^ system for a general 

sweeping 
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fig. 4.11 Example of Cylindrical Sweep Surface 



fig. 4.12 Example of rendering joY extrusion die surface 
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CHAPTER 5 


CONNECTING CURVES AND SURFACES 

5.1 IMTRODUCTION n 

Complicaied die surfaces can be designed by 
conventional method Of supplying control points or data 
points. But the efficiency of the designer will be less. Therefore 
some supporting methods like gluing, blending or intersection 
should be implemented in the software to facilitate the designer. 

Here we attempt a method of connecting B-spline curves 
end to end and connecting B-spline surface patches at the 
boundary. Our method endiures tangency or curvature continuity of 
the connecting curve Or surface patches as required by the 
designer. 

5.B CONNECTING B -SPLINE CURVES 

A complicated ? curve may consist of linear as well as 

k "“2 

curved part. As B-spline curve has everywhere C continuity, Ck 
is the order of the curve 3 an choice of order decides the 
positional, tangency or curvature continuity. For example order S 
joins the control polygon by linear segments, order 3 generates a 
curve inside the control polygon with tangent and order 4 
generates curve with curvature continuity. Most of the die 
surface at best need curvature continuity either on the surface 
or on the boundary curve. 

The following steps are suggested to draw curve 
separately and to connect them end to end. 

1. The maximum order required for any individual curve is chosen 
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as order of the curve. 


2. The whole curve is drawn roughly on paper. 

3. The curve is separated into either linear , curve segmentor a 
combination of both segments. 

4. The control points are input interactively keeping in mind that 
minimum number of control points for each curve should be equal 
or greater than the order of the curve. 

5. If the curves contains both straight line or curved segment then 
there must be two control points at each end of the straight line 
and in addition, one or two control points are added outside 
depending on tangency or curvature continuity required. 

6 The subsequent segments initial control point should coincide 
with the last point of the curve. 

5.3 CONCATEHATI NG CONTROL POINTS AND KNOT VECTOR 

Concatenati ng control is done by combining all the points and 
deleting the redundant points in the array. 

The individual knot vectors for different segment has 
to be concatenated to make it a single curve. The following 
simple algorithm is implemented in our software to combine knot 
vectors into one. 

1. The entire last group of equal values are deleted from the 
first knot vector. 

2. The first element from first group of equal values in the 
second knot vector is deleted. 

3. • n^is added is added to each element in the second knot 
vector, where * n* is the value of the elements removed. 
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4. The knot vector is considered as one. 

Example: 

knot 1 = C O, O, O, 1 » 2, 3, 4, 4 , 4] 
knot a = CO, O, O, 1, 1 , 13 

stepl: knot vector 1 becomes CO, 0, O, 1, 2, 33 
step2; knot vector becomes [0» 0> 1, 1, 13 
stepS: knot vector 2 becomes C4, 4, 5> 5, 53 

step4: The new knot vector is CO, 0,0,1, 2, 3, 4, 4, 5, 5, 53 . 

The above procedure has been implemented in the 
software so that the designer only has to give control points and 
the order for separate segments interactively. 

5-4 SURFACE CONNECTION 
5* 4, 1 Necessary Condition 

Our method connects a boundary of a given B-spline 
surface with the boundary of another surface. The two surface 
and the connecting patch can be unified while keeping tangency 
or curvature continuity everywhere along the joints. 

Suppose two surfaces a and b are given each having at 
least one edge Cnone of them is doubly closedlike a torusD. We 
aim to to construct a patch c which satisfies following 
condi ti ons . 

1. c = cCu,vD is a bi parametric B-spline surface. 

2. c has either plane topology or cylindrical topology. 

3. One of the c^s edges coincides within some tolerance with one 
of a^s edges. The contact curve ca will fix the u-part of the 
parameterization of c: curve ca = cCu,0D. 
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4. Across the boundary between a and c there will be curvature 
CC ^ continuity, although C and C*^continuity can be forced 

whenever necessary. 

5- The speed of varia'tion of "the "tangenL direction or the 
curvature changes continuously along the curve and is kept 
inversely proportional to the local distance between surfaces a 
and b . 

6. The average of this speed variation can be controlled by the 

desi gner . 

7. The condition from 3 to 7 should be also valid for curve cb 

between surface b and c. 

5# 4# 2 Mathematics for Connecting Patches 

The method of generating surface patch for connecting 
two other surface patches proceeds as given belowC133 

1. A net of control points* the B-spline order in u and v 
direction* two knot vectorsC for u and v parameterD are to be 
chosen to describe two surfaces. 

2. c should be closed in the u direction if the relevant edges of 
a and b are closed curves. 

3. The edge cCu*OD will be completely determined by the 
corresponding row of points for c, A number CDZ) of points di Ci 
= DD on the boundary of a should be evaluated and be used as 
interpolation data points, 

4. The algorithm will supply and order geometric continuity that 
ensures lower order and continuity. 
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S. 4-3 Determination of Control Points 

First of all- we determine a unit tangent vector ti and 
radius of curvature p of surface a at d., from which the geometry 
of c can be extrapolated. Fig 5.2CaD. It is sufficient to 
determine the osculating circle of some curve perpendicular to 
the boundary at d.. Point di and two non-degenerate pointsCif the 
points are at least .001mm apartD el and fl on surface a near di. 
fig. 5. 1 are determined approximately with center ml and radius pi 


through di as 

foil ows : 


gx. = ei - 

di 

C5. ID 

hi = fi - 

" di 

C5. ED 

If 

hi 

C5. 3D 

ihi 

C) = 

1^^ 19^1^*" kii^gi + Igi I^C [hi l^-kOhi 



2|gl X hl| 

V. O . 4k J 

1 o 1 = t 

gu 1 |hl I |gl -hi | 

C5. 5D 

1 p . 1 

1. 

2|gl X hl| 

mi = di 4- 

p . Nl 

C5. 6D 




The second term of Eq 5.6 is a vector of length equal to 
radius of curvature p.in the direction of the normal vector Ni 
at dlL103. 

Now a curve cl is tried to fit at point di with start 
direction -gl and with start radius of curvature p^ . One of the 
simplest choice is made for cl is a cubic open B-spline curve 
clCvD with four fold knot vector C £e. order 4D at the beginning Cv = 
03) of its parameter range. The start properties Cat v=0D of 
ciCvD are found in terms of its first 3 control points EK. El, 

Fi : 

clCOD = Du = di C5.7D 

Differentiating the B-spline basis function with respect v and 
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putting V 


O the following values are obtained for a 4th order 


non-rational B-splineCt©] . 

ct = dxdu ciCOD = 3CEi. - DiD C5. 8 ]) 

Cl = d^/ d^u OiCOD = 6 Dl -9El + 3Fl 

= - 6 CEI - DID + 3CFl - ElD C5. 9D 

The radius of curvature at the start of the curve isClOl: 


Rl = 


CiCOD 


ICICOD X CICOD 


3 I El - Dl I 

I Fl - El I si nei 


CS. lOD 


or sin©l = 


3 I El — Dl 


C5. IID 


p ^|Fl - El 


where ©I is the ‘angle between CFl-ElD and CEl - DlD. Further 
constraints are that Dl.-'El and Fl are coplaner with el and fl. 

we have impos'ed the following condition to deter nune 
the point El. 5 


El - Dl 


1 x 4 Rl Cdl - elD 
|dl - el I 


C5. laD 


where Rl is any positive number Cthe choice of Rl will 
be made later in this sectionD. For given Rl, dl,el and flCthus 
given ©iD there still remain several choices for ©i and 
jFl-El| in 3 dimensional field. We have chosen 


|Fl 


Ell 


|El - Dll = 1/4 Rl 


sinei = Rl/p. 

4 t 


C5. 13D 

C5. 14D 


From the above equation point F is determined as follows. 


Referring Fig 5. £CbD 

I . ■ hlD gl 

^ Cg^.gO 

Eh = El + 1/4 Rl cbs©l gl 

IpM 


C5. 15D 
C5. 16D 


Fl 


Eh + 


1/4 "Rl sin©l 


C5. ITD 



5# 4# 4 Choice of conslaint R1 


The detailed development of the curvature of 
CiC star ting with depends largely on Ri. This parameter 

indicates over which range the curvature is influenced by the 
geometry of surface. A simple choice can be the local distance 
between the two surface a and b, where curve Ci intersects them. 

= jCiCmaxD - CiCminZ) | C5.18D 

However » the choice Rl can be varied to have control over the 
connecting surface patch. 

The same procedure described above for surface a must be 
performed on the boundary of b to determine curve Cb at v = 
vamx. The same number of data points are evaluated at the edge of 
edge b. Now we have 6 control points for each dl point on the 
edges. So a surface can be interpolated taking Dx 6 control 
points and of desired degree. 

5* 4. 5 Special cases 

special cases should be treated caref ully^e. g. if di» 
el and fl are col inear then radius of curvature becomes infinite. 
Then the curve Ci will start as a linear interpolation. 

5-5 IMPLEMENTATION 

The above mathematics is implemented in a program 
which joins the v = v-min of second surface to v = vmaLx of first 
surface. Three points are chosen near the boundary to be 
connected very closely to decide the radius of curvature of the 
osculating circle. However, connecting B~spline surfaces can 
hardly be performed exactlyClS]. But the above method 
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approximately meets the condition from 3 to 7 and by deciding the 
number of data points on the boundary of the connecting surface 
and surface patches are found to be reasonably good. 

5.6 USEFULNESS IN DIE SURFACE DESIGN 

As mentioned earlier , the above method can be used as 
toolto increase the efficiency of the surface model by using it 
for connecting arbitrary surfaces obtained by either control 
point manipulation or data fitting and making fillets between 
close or open surfaces or blending two surfaces patches. 

5.7 CONCLUSION 

The quality of connecting patch is checked by only 
gr aphi csC i . e. by line drawings^. In the program the boundary 
fits are not cofined by preset tolerance, nor is, the 
interpolation error evaluated after the fit. Hence,' the 

usefulness for the die design can be realised by a NC product 
model of the shapes. 
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CHAPTER 6 


CONCLUSIONS AND SCOPE OF FUTURE WORK 

6.1 CONCLUSIONS 

In the present work, the main objective was to model 
die sculptured surfaces using various methods available in 

literature for surface modelling, data fitting, sweeping of 

curves to create more complex surface.* In this thesis, an 

attempt has been made to study and implement the above methods to 
describe various die surface geometry. 

First of all it is justified that the choice of 

parametric equation and specifically that of B-“spline is a good 
choice for die surface modelling. 

The second chapter describes a mathematical model for 
data fitting of B~spline surface. In the implementation, less 
than 100 data points with maximum 5th degree splines are 
satisfactorily tried. Theoretically, any number of data points 
to any order of curve could be fitted. However, the dimension 
has been limited to 500 data points and 5th degree curve because 
a curve within 5th degree is sufficient to describe any die 
sculptured surface and secondly Starbase graphics routines allows 
maximum of 5th order splines. The surface fitting technique 
seems to be a very useful method for die surface fitting. The 
interactive design is left for future scope of work. 

Sweep surface technique can really meet the demand for 
variety of open as well '-‘as close die surfaces particularly for 
open forging, die casting, injection moulding and cores for 
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casting. 

Connecting arbitrary B-spline is attempted only to 
increase efficiency of' surface models however » the same is not 
found very useful for the surface modelling of dies. 

The software is presented with menu based program. 
Basically the input is either provided by text file or through 
keyboard. The surface model can be done in wire frame model or 

in solid model. The view can be seen from various angle by 

changing camera position. Some surface model are tried using 

shading and light source. 

The software is incomplete in the sense that it is not 
yet user friendly and interactive design is not implemented for 
all modules. 

6. 2 SCOPE OF FUTURE WORK 

6.2.1 Triiraning - 

Trimming is an important tool for generating complex 
3D surface for a die. 'A trimmed surface is essentially a regular 
tensor product surface "but certain areas are marked as invalid 
or invisible. These ar-eas are usually defined by dense polygons 
within the u,v domain of the patch. Star base graphics has 

trimming of 3D B-spliTie surface routine which can be simply 
defined by control points as described for a surTace. However, 
it does not allow multiple trimming surf ace. Mul ti pie trimming 
could be a future scop^- of work. 

6.2.2 Offset of 3D SctH.ptured Surface 
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Offset of 3D Sculptured Surface 

The offset of the whole 3D surface module is vital for 
die sculptured surfaces for adding tolerances etc. The offset 
technique for 3D surface is still a state of art subJectClQ], 
The offset of the surfaces can be carried out by shifting the 
control polygon in normal direction. But simple shifting control 
points does not necessarily offset 3D surface uniformly. Further 
analysis could be carried out to find out a well defined normal 
of the surface for offset purpose. 

6*2.3 Intersection And Blending 

The intersection of surface patches and consequent 
blending will be very useful for designing various complicated 
surfaces. Particularly^ it will be helpful for the design of 
injection moulding die surface which consist of many intricate 
geometries. This work is being presently carried out by various 
researchersC 203 . 

6*2.4 Machining of parainetric die surface 

To manufacture a die surface defined by parametric 
equation is not very difficult with the advent of powerful CMC 
machines. To machine a 3D surface a minimum 3i/2 axis CNC 
machine is required for milling operation. There are already 
many software packages readily available to generate HC code for 
3D surfaces for example CAM-i , APT> CASPA etc [33. Some possible 
problems involve in real machining conditions are as following 
1* Overcoming the problems of the sharp corners on the surfaces 
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2. Scallop Tolerance - To decide the subsequent course of passes 
of the tool » which should be acceptable between consecutive 
passes. The scallop tolerance decides total work needed for 
surface machining for a finished product. 

3. Gouging- This is a problem arises due to lack of proper 
panning of tool path resulting in under cutting in some part. 

4. Adaptive path planning- Which considers local optimum paths 
which are assembled to a continuous global path. 

5. Tool and work piece interference: In complex 3D surface 
machining tool might interfere with the work piece unless path 
checked clearly. 

Taking into account of the above problems the die 
surface machining by CNC machines for either to produce the die 
itself or to use EDM could be further point of research for the 
development of integrated CAD/CAM of dies. 
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APPENDIX 


USER' S MANUAL 



ESC This key is 10 be pressed bo reburn bo previous menu. 

curser keys are bo be used for moving inside bhe menu. 

SELECT: Ib selecbs various sub-menu as lisbed above. Any ibem can 
be selecbed by pressing [ reburn] key. Afber selecbion, bhe 
necessary value is asked on bhe righb hand corner, for example if 
inpub mode and file are selecbed ib asks for file name. If 
inberacbive is selecbed ib asks for ail values required bo make a 
file. A bexb file is required in bhe. form as asked by bhe 
inberacbive mode. 

DISPLAY: Ib displays bhe final graphics. To reach bhis menu ESC 
key is bo be pressed. 

EXIT: When pressed ib exibs from bhe main menu. To exib from a 
bexb file '^C is bo be pressed bwice. 
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